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We review some aspects of recent work concerning double scaling limits 
of singularly perturbed hermitian random matrix models and their connec¬ 
tion to Painleve equations. We present new results showing how a Painleve 
III hierarchy recently proposed by the author can be connected to the 
Lenard recursion formula used to construct the Painleve I and II hierar¬ 
chies. 
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1. Introduction 

The semi-classical one hermitian matrix model is defined as the measure 

1 g-nTrI/(M)^^ (1) 

on the set of n x n hermitian matrices whose spectrum is a subset 

of intervals J. Here the potential V is such that V' is a rational function. 
The angular degrees of freedom of M may be integrated out of this model 
to give a j.p.d.f for the eigenvalues Xi of M, 

1 ” 

( 2 ) 

where w{x) := ^ A is the Vandemonde determinant and XJ b the 

indicator function on J. Such models are a very general class of models for 
which the method of orthogonal polynomials can be used to give a solution. 
The method of orthogonal polynomials expresses the eigenvalue /c-point cor¬ 
relation functions in terms of the correlation kernel, 

Kn{x,y) = (p^(x)pn-l{y) - Pn{y)Pn-lix)) , (3) 

X y 


( 1 ) 
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where Pj, j = 0,1,... are a family of monic polynomials of degree j charac¬ 
terised by the relations 

J Pj{x)pmix)w{x)dx = hj5jm- (4) 

The limiting mean eigenvalue density is given by 

p{x) := lim -Kn{x,x), (5) 

n—¥oo 71 

and describes the macroscopic behaviour of the eigenvalues for large n. 

Such models have been studied at finite n in [T] and their relation to 
integrable systems fully described. What is less known is the types of critical 
behaviour in such a model. In the case that V is polynomial and J = M, a 
number of distinct critical points have been identified and studied over the 
last two decades. These have been classified as; 

• Edge: The spectral density acquires extra zeros at an edge a of its 
support. The usual behaviour of p near a is p{x) = 0{\x — a\"^) 
however when extra zeros are present the possible behaviours of p 
are p{x) = 0{\x — a\^~) with k € N®. The limiting kernel in the 
neighbourhood of such points is constructed in terms of solutions to 
the 2A:-th Painleve I equation. 

• Interior: The spectral density acquires new zeros at some point a in the 
interior of its support. The behaviour of p near a is p{x) = 0{\x — a\^^) 
with k € N. The limiting kernel in the neighbourhood of such points 
is constructed from solutions to the fe-th Painleve II equation. 

• Exterior: The spectral density acquires a new cut in its support. This 
transition is known as “birth of a cut” [8]. Here the limiting kernel is 
constructed from Hermite polynomials and the local behaviour in the 
new cut mimics a GUE matrix model. 

In the more general case of the semi-classical model no such classification 
exists, however a number of special cases have been investigated in the 
recent literature. 

• The effect of logarithmic singularities in the potential have been inves¬ 
tigated in a number of works. In [7] the effect of a singularity in the 
bulk results in a kernel constructed with Bessel functions. In [9] the 
situation of a logarithmic singularity coinciding with a interior critical 
point was found to lead to kernels containing solutions to the general 
Painleve II equation. Einally logarithmic singularities at the edge of 
the spectrum results in general Painleve XXXIV equations. 
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• The addition of a hard edge also results in new behaviour. It has long 
been known that the kernel near the hard edge can be constructed 
in terms of Bessel functions. More recent work m has considered 
the case of a hard edge meeting a soft edge, with the resulting kernel 
constructed in terms of Painleve XXXIV transcendents. This was 
further extended in m to a hard edge meeting a edge critical point. 
It was shown that the associated Painleve transcendents satisfy the 
k-th member of the Painleve XXXIV hierarchy. 

• Finally, very recently, the behaviour of eigenvalues near poles in the 
potential have been studied. The case of a simple pole both in the bulk 
and at the hard edge have been investigated in lasiEiis] and it was 
shown that the kernel can be constructed using solutions of Painleve 
III. The case of higher order poles at a hard edge was studied in [6] 
by the author and collaborators. The kernel in the neighbourhood of 
the pole was constructed using solutions of a member of a Painleve 
III hierarchy. 

In this short note we report on some new aspects of the work in [6]. In 
particular we give a relation between the Painleve III hierarchy defined in 
[6] and the Lenard recursion relations that are ubiquitous in the Painleve I 
and II hierarchies. 


2. A Painleve III Hierarchy 

In [6] the fc-th member of the Painleve III hierarchy was defined as the 
system of k ODEs {p = 1,... ,k), 


p 

^ ^ if-k—p+q+l^k—q ~ i^k—p+q^k—q) + ^^k—p+q^k—q ~ ^'^^k—p+q^k—q) — Lp) 
q=0 

( 6 ) 

for k unknown functions ii = ... = ^k{s), with = 0 and 

io{s) = |. The Tp’s are constants that act as times. The quantity u = u{s) 
is defined by, 

u{s) = -^ ((4)" - 3(4)^ + ^o) • (7) 


Example 2.1 For k = 1 we have the equation 

= 4('g)^ _ 4('S) _ _ To n 

^ ll{s) S S ^l{s) s 


(8) 


which we identify as a special case of the Painleve III equation. 
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Fig. 1. The jump contour F for the model RH problem for dt. Contours are labelled 
by S and sectors by fl. 


Example 2.2 If /c = 2 we have a system of two ODEs; 

n _ To l'i{s) e^{s) l2{s) 

2ii{s)i2{s) l2{sY i2{sY £i(s)£2(s) 4('S) ^2('S) 2£i(s) 


£i(sff2(sf 






STo 


Us) 


UU Us) 


- T2 


2Us)^Us) 

Us) 


- 2li{s)l'i{s) + sl'i{s) + 2 l 2 {s)h{s). 


( 10 ) 


3. A Riemann-Hilbert Problem for the Painleve III Hierarchy 

In [6] it was shown that a solution to the fc-th Painleve III equation may 
be extracted from the following RH problem. 

(a) <I> : C \ S —> analytic. See Figured] 

(b) <I> has the jump relations ^+(-s) = ^-{z)ji for z G Sj 



(c) As z —)• oo, <1> has the asymptotic behaviour 

1 / ri;(s) u(s) 


^z) = 


1 0 

v{s) 1 


I+- 


Z \h{s) -w{s) 


+ 0(zU 


1/2^ 


X ' <"3, (12) 


where + ’^^i) v, h and w are functions of s. 















lenard printed on February 18, 2015 


5 


(d) As 2 ; —>■ 0, there exists a matrix <ho(s), independent of z, such that 
has the asymptotic behaviour 


<^{z) = Ms){I + 

for 2 : G ilj, where Hi, H 2 , are given by 

Hi = 1, H2= fj, Hs = 



(13) 


(14) 


We then have as a corollary of [B] Theorem 1, that 

Theorem 3.1 Let a > —1, and let 4>(2;;s) he the unique solution of the 
above model RH problem for s > 0. Then, the limit 


vis) 



lim s$(z, 

.z—^oo 


21 


(15) 


is a solution of the k-th member of the Painleve III hierarchy 


Remark 3.2 The proof of this theorem identifies y with ii in the Painleve 
III hierarchy. Furthermore it also demonstrates that the Lax pair for <I> 
takes the form, 


A( 2 :, s ) = a{z, s)a 3 + b{z, s)(T+ + c{z, s)a-, (16) 

B{z-, s) = {z — u{s))a-+ a+. (17) 

where o, b and c are related by, 

a{z,s) =-^dsb{z,s), (18) 

ciz, s) = {z- u)b{z, s) - ^d‘^b{z, s), (19) 

dsc{z, s) = 1 + 2{z — u{s))a{z, s). (20) 

Substituting (fT8|l and (fT^ into (12(1 yields 

zdsh{z, •s) = ^ {dlh{z, s) + ‘iu{s)dsb{z, s) + 2u'{s)b{z, s)) + (21) 


We may compute b{z,s) by substituting b{z,s) = Ylj=o^k-jis)i4:zy 

into the above equation to get, 

^'j+iis) = + 4:u{s)i'j{s) + 2u'{s)ij{s), 


( 22 ) 
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4. Integration of Lenard-type Recursion Relations 

The recursion relation (|22p is the Lenard recursion relation appearing 
in the Painleve I and II hierarchies. In those cases the initial condition 
is ^0 = ^ whereas here we have £o = 'S/2. Let us consider the general 
case where and are known functions. The fact that is known 
implies u{s) satisfies a integro-differential equation of order 3A: + 1. The 
following lemma gives k + 1 constants of motion, i.e. functions of u{s) and 
its derivatives which are constant in s. This allows the equation for u{s) to 
be reduced to an ODE of order 2k and we will see that these constants of 
motion are precisely the ODEs appearing in ([6|). 


Theorem 4.1 Let ij he the integro-differential polynomials in u generated 
by the Lenard recursion relation (I22p together with an initial condition for 
(.Q. Furthermore let (-k+i have a given form. The integro-differential 
equation corresponding to dk+i I'he following constants of motion; 


T~p — ^k+ffk—p T 'y ^ {f-k—q^k—p+q+i ^k—p+q,k—q) j 0 ^ p ^ fc (^3) 
q=0 

if ^'k+i ~ 


p-i 

(Tp = (-O^p + y ^ {fdp—X—q,q (-p—X—q^q+l) i 0 < p < /c (24) 

q=0 

if iQ = 0. In the above expressions we have introduced, 

■= AA” - 3^4 + (25) 


Proof. We begin with the following identity, 

+ ^n(-m+l ~ ^n,m (26) 

This identity can be established by the following argument. 


AAlY = ^mC + + 2n'^™4 (27) 

= £AY + " UC+1 " " C) (28) 

= imC + inC + Huimin)' " 44+1 (29) 

£m£n+l + 44+1 = ( 30 ) 


where the first line is the Lenard-type recursion relation multiplied by 4) 
the second line is obtained by grouping terms and applying the recursion 
relation for im- The final line is obtained from the elementary identity, 
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Now note that integrating (l26li by parts and letting n n — 1 gives, 

• (32) 

Geometrically this equation says that the quantity Imt-'n only picks up total 
derivatives if we move along the anti-diagonals of the lattice points labelled 
by (m, n) for 0 < m,n < k + 1. We can now see how the constants of motion 
arise. If, due to the boundary conditions on Iq and Ik+i^ the quantity 
is a total derivative on the border of the lattice (i.e. when m or n are 0 or 
A; -|- 1) then we can produce a total derivative equal to zero by using 
to move across the lattice from one border to another. Explicitly, moving a 
distance r along an anti-diagonal gives. 


/ f' 


i' + 

^n.—r ' 




r—1 


^ ^ ^n—g—l,m+g ^n—q—1^7n-\-q+l 

q=0 


(33) 


If ^'k+i — 0 then let m = A: — p, n = A: -|- 1 and r = p +1. This gives. 


d 

ds 


p 

^k+l^k—p T ^ ^ (life 

q=0 


—p+q,k—q 


^k—q(-k—p+q+l) 


= 0 


(34) 


which integrating gives the first part of the theorem. If = 0 then let 
m = 0 and r = n, this gives, 


ds 


p-i 

^O^p ^ ^ (S^p—l—q,q ^p—l — q^q+l) 

q=0 


= 0 , 


(35) 


which integrating gives the second part of the theorem. 


□ 


Remark 4.2 Setting ik+i = 0 and = s/2 in the Lenard recursion relation 
implies the Tp are constant and we recover (0. 

Remark 4.3 Consider the standard Lenard differential polynomials ob¬ 
tained with the boundary condition Iq = ^ and setting all integration con¬ 
stants to zero. From the above theorem we find that ap are constants, which 
by the definition of the standard Lenard differential polynomials must be 
zero. We therefore have, after some rearranging, 

p-2 

^p — ^ ^ (Up—1— ('p—l—q('q+l') T fIo,p—!■ 

q=0 


(36) 
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The right hand side of the above equation only contains In with n < p while 
the left hand side only contains ip. Hence we can use these equations to 
recursively determine ip, which are exactly the Lenard differential polyno¬ 
mials. Note that this shows each Lenard differential polynomial is indeed a 
differential polynomial, a fact that is not obvious from their usual definition 
(I22D. 
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